Full counting statistics for the Kondo dot. 
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The generating function for the cumulants of charge current distribution is calculated for two 
generalised Majorana resonant level models: the Kondo dot at the Toulouse point and the resonant 
level embedded in a Luttinger liquid with the interaction parameter g = 1/2. We find that the low- 
temperature non-equilibrium transport in the Kondo case occurs via tunnelling of physical electrons 
as well as by coherent transmission of electron pairs. We calculate the third cumulant ('skewness') 
explicitly and analyse it for different couplings, temperatures, and magnetic fields. For the g — 1/2 
set-up the statistics simplifies and is given by a modified version of the Levitov-Lesovik formula. 
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Since Schottky's realisation that the shot noise in 
a conductor contains invaluable information about the 
physical properties of the charge carriers, the question 
about the noise spectra of different circuits became as 
important as the knowledge of their current-voltage char- 
acteristics 0. The noise constitutes the second moment 
of the current distribution function (which is the prob- 
ability of measuring a given value of the current) and 
is supposed to contain information about the charge of 
current carrying excitations at weak transmission (reflec- 
tion). That is indeed the case for S — S and S — N junc- 
tions LI but not been proven for ageneric interacting 
model. As has been pointed out in |fj, the third cumu- 
lant also contains valuable information about the charge 
of the current carriers. Therefore it is natural to investi- 
gate the full current distribution function. This was an 
academic question for a very long time as even the mea- 
surement of the second cumulant remained on the frontier 
of experimental physics. Only after the work by Reulet 
and coworkers the measurement of the third cumulant 
became possible Q . Inspired by this remarkable achieve- 
ment, the full current distribution function (more often 
referred to as 'full counting statistics' or FCS) has been 
theoretically analysed in recent year for a wide range of 
systems. 

In their seminal work Q , Levitov and Lesovik derived 
the exact formula for the FCS for the electron tunneling 
set-up (single channel): 

hx Xo (X;V;{T(cj)})=T jf° ^l n {l + T( W ) (1) 

x [n L (l - n R ){e lX - 1) + n R (l - n L )( e -* x - 1)] } , 

where T is the waiting time, A is the measuring field 
(we will explain this notation in more detail shortly), 
n H (i)(o;) = tif(u) ± V/2) are the electron filling factors 
in the (right and left) leads, V is the bias voltage, and 
T{u>) is the single electron transmission coefficient. The 
knowledge of T(u>) thus fully defines the FCS for non- 
interacting systems. 

While the Levitov-Lesovik approach can be relatively 



easily generalised to various multi-terminal (multi- 
channel) set-ups, it is notoriously difficult to include 
electron-electron interactions. Up to now most works 
in this direction relied upon variousperturbative expan- 
sions, in the tunnelling amplitude [a, or in the inter- 
action strength (for a recent review see ^3), as WC U 
as calculations at zero temperature [TTI |. Certainly no 
paradigm for an interacting FCS em erge d as of yet. No- 
table exceptions are the contribution \12\ as well as recent 
works by Andreev and Mishchcnko and Kindermann and 
Trauzcttcl [Hill, where the exact FCS was calculated 
for the (single-channel) Coulomb Blockade (CB) set-up 
of Matveev and Furusaki 0, ^| . We shall establish the 
precise connection of these works to our results. 

The purpose of this Letter is to contribute to our un- 
derstanding of interacting FCSs by means of obtaining 
the exact FCS for two particular experimentally relevant 
set-ups: the Kondo dot and the resonant tunneling (RT) 
between two g = 1/2 Luttinger liquids (LL). 

We start with a brief description of the method. Our 
goal is the calculation of the generating function x(A) = 
'^2e' lqX Pg for the probabilities P q of q electrons being 
transmitted through the system over time T. We first 
define operators Tj transferring one electron through the 
system in the direction of the current (i = R) and in the 
reversed direction (i = L). The electron counting opera- 
tor on the Keldysh contour C can then be written down 
in its canonical form as 

T x = e iX W 2 T R + e- iX W 2 T L , 

where the measuringing field X(t) is explicitly time de- 
pendent, A(i) = X9(t)6(T — t) on the forward path and 
A(i) = —X9(t)9(T — t) on the backward path. Accord- 
ing to |a, [7| , the generating function is then given by the 
expectation value 

/ -ijT x (t)dt\ 

xW = {T c e c \ . ( 2 ) 

In order to calculate x(X) we define a more general func- 
tional x[X-(t), X+(t)] formally given by the same Eq. @ 
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but where X(t) is now understood to be an arbitrary func- 
tion on the Keldysh contour, ± referring to two different 
functions on the time and anti-time ordered halves of the 
contour. Next we assume that both X±(t) change slowly 
in time. Then, neglecting switching effects, one obtains 
at large T 

X [X-(t),X+(t)]=e X p{ -i [ U[X^(t),X + (t)]dt\ , (3) 



where Z-/(A_, A+) is the adiabatic potential. Once the adi- 
abatic potential is computed, the statistics is recovered 
from ln%(A) = —iTU(X, —A). As X± are external param- 
eters, performing the derivative of both J3J and with 
respect to, say, A_, with help of the Feynman-Hellmann 
theorem |f7j ] we immediately obtain 



d 
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U(X-,X+) = 



where we use notation 



X(A-,A + ) 



T c 



9A_ l x 



-ijT x (t)dt 

A(t)e c 



This is somewhat more complicated than the usual 
Hamiltonian formalism for a quasi-stationary situation, 
we'll give further technical details in the long version [Tsj : 
in particular, we have verified that Eq. comes out cor- 
rectly in the non-interacting case. 

In order to study the FCS for the Kondo dot we use 
the bosonization and rcfermionization approach, origi- 
nally applied to this problem by Emery and Kivelson [l9| 
(see also |2(]]) and refined by Schiller and Hershfield (SH), 
see [2lJ. The starting point is the two-channel Kondo 
Hamiltonian (we set h = Vf = e = ks = 1 throughout), 



H = H Q + H, + H. 
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where, with ^} a ,a being the electron field operators in the 
R,L channels, 



H 
Hj 

H v 
Hm 



i H / dx^ aa {x)d x ij: aa {x) . 

a=R,£ <r=T,J. J 

J " /3s a/3 T!/ ' 



a,p=R,L v=x,y,z 



{V/ 2) W ' dx (iPl^La ~ , 

-ia Bgi HT z = -At 2 . (4) 



Here t v=x ' v ' z are the Pauli matrices for the impurity spin 
and (a, [3 — R 1 L; a =| , |; a v aa , are the components of the 
fth Pauli matrix) 

*a/J = E^L(0)^^(0), 



are the electron spin densities in (or across) the leads, 
biased by a finite V. The last term in Eq. J3J stands 
for the magnetic field, A = ^,B9iH . Following SH, we 
assume = Jf = jf, J 2 iL = jf R = J z and J 2 ii? = 
J Z L = 0. The only transport process then allowed is the 
spin-flip tunnelling, so that we obtain for the T\ operator 



J 



RL 



+ p iA(t)/2 ,t 



r 1 e 



After bosonization, Emery-Kivelson rotation, and 
refermionization (see details in plj l and going over to 
the Toulouse point J z = 2n, which is the only approxi- 
mation we make, one obtains 

H' = Hq — i(J_ b£f + J + ar)f) — iAab + T\ , (5) 

where the counting term is given by 



T x = -iJ ± b K cos(A/2) - rj sin(A/2)] 



(6) 



with J± = (Jf L ± Jf fl )/V27ra^ J± = Jf i / v / 2iW (a 
is the lattice constant of the underlying lattice model) 
and a and b being local Majorana operators originating 
from the impurity spin. The fields r/f and in the spin- 
flavour sector are equilibrium Majorana fields, whereas rj 
and £ in the the charge-flavour sector are biased by V, 



H' a = i I dx r lf {x)d x r lf {x)+i f {x)d x £,f{x) (7) 



+ r,(x)d x r,(x) + £(x)d x S(x) + VZ(x)t)(x) 



Using Eqs. ©-0 on e can straightforwardly evaluate the 
adiabatic potential W(A_, A+) as the problem has become 
quadratic in the Majorana fields. 

Skipping details of the calculation, we report the re- 
sulting exact formula for the FCS of the Kondo dot, 
which is the main result of this paper: 
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lnx(A) = T f ^ ln{l + 7\(u;) [n L (l - n R )(e 2lX - 1) + n R (l - n L )( e - 2lX - 1)] (8) 
o 

+T 2 {uo) [[n F (l - n R ) + n L (l - n F )](e lX - 1) + [n F (l - n L ) + n R (l - n F )( e - lX - 1)]] } , 

where now the filling factors are n R ^ = n F (u> ± V) [n F (uj) being the conventional Fermi function], not to be confused 
with notation in Eq. |T|l. The 'transmission coefficients' are 



„ rj(^ 2 + r 2 + ) 2r ± r_(^ 2 + r 2 + ) + 2A 2 r ± ^ 



[c^ 2 - A 2 - r+(r ± + r_)] 2 + c^ 2 (r + + r_ + r ± ) 2 [ w 2 - a 2 - r+(r ± + r_)] 2 + ^ 2 (r+ + r_ + r_ 



\2 ' 



where T 4 = J 2 /2 (i = ±,_L). 

For small voltages and zero temperature the generating 
function turns out to be quite simple. 



X (A) = [1 + T e (e* x -1)] 



(9) 



where TV = TV/ir is the number of the incoming parti- 
cles during the time interval T and T e = ^/Xi(0) is the 
effective transmission coefficient. We analysed explicitly 
the behaviour of the system around the Toulouse point. 
The trivialisation © turns out to be robust against de- 
parture from this special point in parameter space |l8j| . 
While (|9j| is the standard Lcvitov-Lesovik result for spin- 
ful systems, the physical content of Eq. (jHl is more inter- 
esting, since it cannot be reduced to binomial statistics 
as in Eq. 10 at finite T. One can see that the charge 
current is carried by two different quasi-particles with 
charges qi = 2e and q 2 = e. We identify the correspond- 
ing transmission coefficients as T\ and T 2 , respectively. It 
was realised by SH that at least in the case of finite mag- 
netic field A ^ it can become energetically favourable 
to tunnel electr on p airs through the impurity rather than 
single electrons [21| . The presence of the term containing 
2 A in © can be interpreted as a signature of this effect. 
The full transport coefficient To as calculated by SH turns 
out to be a composite one and it is recovered from T\^ 2 
through a very simple relation: T = Ti+T 2 /2. From the 
point of view of the Kondo physics, the case when T 2 = 
and the statistics reduces to a modified Levitov-Lesovik 
formula, x(A) = Xo /2 (2A, 2V, {Ti(cj)» (binomial statis- 
tics at T = 0), is the symmetric model in zero field (the 
other, unphysical, case of T 2 = is when J± =0). We 
have evaluated the first and the second cumulant of the 
Kondo FCS Eq. (jHlwhich are the same as calculated by 
SH at all V and T jig] . We shall not reproduce these two 
cumulants here and concentrate instead on new results. 

The full analytic expression for the third cumulant ex- 
ists but is too lengthy to be given here. We shall rather 
investigate various limits and use numerics for the gen- 
eral case. So, at T = we obtain: 

(Sq 3 ) = T — [T 2 + 8Tj - 3(T 2 + 2Ti)(T 2 + 471) 
Ja 2tt 

+ 2(T 2 + 2T0 3 ] . (10) 



In the zero magnetic field it yields the following limiting 
forms: 

<^ 3 )^o « TGo ^-^J^V , (11) 
(8q 3 ) v^oo « T7rG r ± , 

where Go = l/(27r) is the conductance quantum. At low 
voltages the cumulant is negative for T_ < r_i_. Gener- 
ally, under these conditions the n-th cumulant appears 
to possess n — 2 zeroes as a function of V, according to 
numerics. The saturation value in the limit V — ► 00 is 
independent of the coupling in the spin-flavour channel 
because the fluctuations in the biased conducting charge- 
flavour channel are much more pronounced than those 
in the spin-flavour channel, which experiences only rela- 
tively weak equilibrium fluctuations. 

In the opposite case of near equilibrium all odd cumu- 
lants (Sq 2n+1 ) are identically zero, which can readily be 
seen from Eq. (JSJ by substituting tirl = n F into it. In 
the limit of low temperatures T — » we recover the con- 
ventional Johnson- Nyquist noise power |2l]|. Moreover, it 
can be shown that the leading behaviour in temperature 
of every even order cumulant in this situation is linear, 
e. g. for (Sq 4 ) we obtain 



(Sq 4 ) « T4G T 



r ± r+(r_r H 



[A 2 + r+(r ± + r_)r 



(12) 



For the general situation of arbitrary parameters, the 
cumulants can be calculated numerically. The asymp- 
totic value of the third cumulant at high voltages, simi- 
larly to the findings of 0], does not depend on temper- 
ature and is given by the result l|llf) . see Fig. ^ 111 the 
opposite limit of small V, (Sq 3 ) can be negative. Suffi- 
ciently large coupling T_ or magnetic field suppress this 
effect though [T^ . 

According to the result of Ref . 0] , as long as the dis- 
tribution is binomial, (Sq 3 ) / (Sq) = (e*) 2 , where e* is the 
effective charge of the current carriers. This quantity is 
to be preferred to the Schottky formula because of its 
weak temperature dependence. Indeed we find numeri- 
cally that the ratio (Sq 3 )/(Sq) in the present problem is 
weakly temperature dependent (it is flat and levels off to 
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(unphysical case). Consequently, the FCS is given by a 
modification of the Levitov-Lesovik formula: 



FIG. 1: The voltage dependence of the third cumulant for 
different temperatures and zero magnetic field (A = 0) for 
r_/r_L = (main graph), and for r_/r± = 0.9 (inset). 



1) in comparison to (Sq 2 ) / (6q) , which is indeed in accor- 
dance with 0. 

We now briefly turn to the g = 1/2 RT set-up. This 
set-up has caused much interest recently, see Ref. [2^ 
and references therein. The Hamiltonian now is 

H = H + + d^ R + H.c) + AdU + H c , (13) 

where Ho stands for two biased LLs at g = 1/2, d is the 
electron operator on the dot, 7 is the tunneling amplitude 
and He is an electrostatic interaction we do not write 
explicitly here (see H^)- Introducing A as standard, and 
carrying out the bosonization-refermionisation analysis, 
we find the same set of equations as for the Kondo dot, 
Eq. © and Eq.©, but with A -> A/2 and J± = 2-f, 
J± — 0, when the Kondo statistics simplifies to binomial 



Xi/2(A)=Xo /2 (A;2V;{T A H}), 



(14) 



with the effective transmission coefficient 7a (w) = 
47 4 w 2 /[47 4 w 2 + (lu 2 - A 2 ) 2 } of the RT set-up in the sym- 
metric case |22j (the contact asymmetry is unimportant). 
All the cumulants are thus obtainable from those of the 
non-interacting statistics Eq. QJ. 

The A = RT set-up is equivalent to the model of 
direct tunneling between two g = 2 LLs Jj&j. The lat- 
ter model is connected by the strong to weak coupling 
(1/g — > g) duality argument to the g = 1/2 Kane and 
Fisher model which is, in turn, equivalent to the 
CB set-up studied by KT. Therefore their FCS must be 
related to our Eq. I|14|) at A = by means of the trans- 
formation: To — > 1 — To and V — ► V/2. Indeed after some 
algebraic manipulation with KT's Eq. (12), we find that 
the FCS for the CB set-up can be re-written as: 

Xcb{X)=^ /2 {-X;V;{1-T {u)}). 

To summarise, we derived the generating function for 
the charge transfer statistics for the Kondo dot in the 
Toulouse limit and analysed the third cumulant in de- 
tail. At low temperatures the transport is accomplished 
by electrons as well as electron pairs in the generic case 
whereas at T = the conventional binomial statistics is 
restored. 
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